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Abstract 

High temperature reduction of the SU(2) Higgs model is realised 
by partially integrating its partition function. Various approximate 
forms of the effective theory resulting from the integration over non- 
static fields and the static electric potential are analysed. Also non- 
polynomial and non-local terms are allowed. Consistency of the per- 
turbative solution is ensured by new types of induced counterterms. 
Perturbative phase transition characteristics are presented in the Higgs 
mass range 30-120 GeV, and compared to results of other perturbative 
approaches. 
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1 Introduction 



Physical systems develop at finite temperature multiple mass-scales. Decou- 
pling theorems Q ensure that effective representations can be found for the 
lighter degrees of freedom at scales asymptotically separated from the heav- 
ier scales. The effective couplings can be related to the parameters of the 
original system by matching values of the same quantities calculated in the 
effective model and the original model simultanously || [|, Ell- 
in weakly coupled gauge theories (the class where the electroweak the- 
ory belongs to) matching can be performed perturbatively. The simplest is 
to match the renormalised Green's functions computed in the two theories 
separately. This tacitly assumes the renormalisability of the effective model. 
Though at large distances local operators of lower scaling dimensions domi- 
nate indeed, it is important to assess quantitatively the rate of suppression 
of higher dimensional operators. 

In finite temperature theories integration over non-static fields yields for 
static n-point functions with typical momenta 0(p), expressions analytic in 
p 2 /M% eavy , therefore the occuring non-localities can always be expanded in a 
kind of gradient expansion. We shall demonstrate, that in Higgs systems to 
0(g A , Xg 2 , A 2 ) accuracy, the non-local effects in the effective Higgs-potential 
arising from non-static fields can be represented exactly by correcting the 
coefficient of the term quadratic in the Higgs-field. 

Similar conclusions are arrived at when the contribution of non-static 
modes to higher dimensional (non-derivative) operators is investigated. 

The situation changes when a heavy static field is integrated out. The 
contribution to the polarisation functions from this degree of freedom is not 
analytic neither in p 2 /M% e nor in & 2 /Ml eavy . For the momentum depen- 
dence no gradient expansion can be proven to exist. Though the contribu- 
tion to the effective Higgs-potential can be represented approximately by an 
appropriately corrected local gauge-Higgs theory, it does not automatically 
correspond to a systematic weak coupling expansion. Assuming that the 
important range of variation of the Higgs field is ~ gT, one proceeds with 
the expansion, and the validity of this assumption is checked a posteriori. In 
case of strong first order transitions this assumption is certainly incorrect, 
but seems to be justified for the Standard Model with Higgs masses around 
80GeV. 

Non-analytic dependence on field variables is reflected in the appearence 
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of non-polynomial local terms in the reduced effective theory. Consistent 
treatment of such pieces is an interesting result of our paper (see the discus- 
sion of the U(l) case in JTDJ). 



The method of integration over specific (heavy static and non-static) 
classes of fields in the path integral (Partial Path Integration, PPI), instead 
of matching, provides us with a direct method of inducing all non-local and 
non-polynomial corrections automatically. Though the calculations are much 
more painful, than those involved in matching, the value of PPI in our eyes 
is just the possibility of assessing the range of applicability of the technically 
simpler method. 

The explicit expression of the non-local and/or non-polynomial parts of 
the action and the impact of these terms on the effective finite temperature 
Higgs-potential will be investigated in the present paper for the SU(2) gauge- 
Higgs model. This model is of central physical importance in investigating 
the nature of the cosmological electroweak phase transition [||]. Earlier, PPI 
has been applied to the O(N) model by one of us ||. 

A similar, but even more ambitious program is being pursued by Mack 
and collaborators |7|]. They attempt at performing PPI with the goal of 
deriving in a unique procedure the perfect lattice action of the coarse grained 
effective model for the light degrees of freedom. 

Our calculation stays within the continuum perturbative framework. PPI 
will be performed with two-loop accuracy for two reasons: 

i) The realisation of PPI on Gaussian level does not produce non-local 
terms, though non-polynomial local terms already do appear. The consistent 
cut-off independence of the perturbation theory with non-polynomial terms 
can be tested first in calculations which involve also 2-loop diagrams. 

ii) It has been demonstrated that the quantitative non-perturbative char- 
acterisation of the effective SU(2) gauge-Higgs model is sensitive to two- loop 
corrections of the cut-off (lattice spacing) dependence of its effective bare 
couplings H |S|. We would like to investigate possible variations of the rela- 
tion of the effective couplings to the parameters of the original theory when 
the applied regularisation and renormalisation schemes change. 

The integration over the non-static and of the heavy static fields will not 
be separated, but is performed in parallel with help of the thermal countert- 
erm technique already applied in the PPI of the U(l) Higgs model JTU . 



The paper is organised in the following way. In section 2 the basic defini- 
tions are given and the operational steps involved in the 2-loop PPI are listed. 
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In Section 3 we are going to discuss the contributions from fully non-static 
fluctuations. In Section 4 the contribution of the diagrams involving also 
heavy static propagators is discussed. In both sections particular attention 
is payed to the analysis of the effective non-local interactions. In Section 
5 the effective Higgs potential is calculated from the 2-loop level effective 
non-local and non-polynomial 3-d (NNA) action with 2-loop accuracy. Also 
the local polynomial (LPA) and local non-polynomial (LNA) approximations 
are worked out. The quantitative perturbative characterisation of the phase 
transition and its comparison with other approaches will be presented in Sec- 
tion 6. Our conclusions are summarised in Section 7. In order to make the 
paper better readable, the details of the computations discussed in Sections 
3 to 5 are relegated to various Appendices. 

2 Steps of perturbative PPI 

The model under consideration is the SU(2) gauge-Higgs model with one 
complex Higgs doublet: 

£ = iKnKn + §(V™$ ) t (V m $) + |m 2 $t$ + A($t$)2 + £cT? 

Fran = ^mA® — ^n^m + 9 e ab C A b m A c n , 

V m = d m - igA a m r a . (f ) 

The integrations over the non-static Matsubara modes and the static electric 
component of the vector potential will be performed with two-loop accuracy. 
In the case of the static electric potential resummed perturbation theory is 
applied, simply realised by adding the mass term (Tm 2 D /2) J d 3 x (Ag(x)) 2 
to the free action and compensating for it in the interaction piece. The 
resummation discriminates between the static and non-static Aq components, 
therefore in the Lagrangean density the replacement 

A (x,r) -> A (x) + a (x,r) (2) 

should be done. (With lower case letters we always refer to non-static fields.) 

In the first step we are going to calculate the local (potential) part of the 
reduced Lagrangean. For this a constant background is introduced into the 
Higgs- multip let: 

$(x,r) = $ o + 0(x,r), 
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0(x,r) 



(3) 



where $ was chosen to be the constant counterpart of £4. The dependence on 
the other static components is reconstructed by requiring the 0(4) symmetry 
of the potential energy. 

In the second step the kinetic piece of the reduced Lagrangean is investi- 
gated. In order to fix the coefficient of the conventional kinetic piece to 1/2, 
one has to rescale the fields. The wave function renormalisation constants 
should be calculated with 0(g) accuracy. The derivative part of the scalar 
action can be extended upon the requirement of spatial gauge invariance to 
include the magnetic vector-scalar interaction into itself. 

Higher derivative contributions to the effective action can be summarized 
under the common name non-local interactions. They appear first in the 
two-loop calculation of the polarisation functions of the different fields. 

The exponent of the functional integration has the following explicit ex- 
pression: 



£(cO + £(2) = 



-m 2 <$>l + 
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K K 



<{K) 



K 2 + m 2 + 



A 



-<t> 



+ 2 U 



-K) 



K + m 



2 + ~^l 



UK) + ci(K)K 2 c a (K). (4) 



where c denotes the ghost fields, and the notation £ 4 = £ H has been intro- 
duced. Landau gauge fixing and the Faddeev-Popov ghost terms are included. 
The above expression implies the following propagators for the different fields 
to be integrated out: 



(a a m (P)a b n (Q)) 



Jab 

P' 2 + m'i 



p p 

X 1 m 1 n 



S(P + Q), 



(A« (P)A b (Q)) 



0~ab 



P +rn A0 



^S(P + Q), 
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2 , 9 ^2 

m D + — $ 



o- 
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(UP)UQ)) = p>l m >J ( P + m 2 H = m 2 + ±& Q , 

(UP)UQ)) = p^T*( P + Ql ml = m* + ^l, 

(cl(P)c b (Q)) = - 5 ^6(P-Q). (5) 

(5 denotes the finite temperature generalisation of Dirac's delta- function.) 

The interaction part of the Lagrangean density consists of two parts. In 
the first set all non-quadratic vertices are collected: 

C x = igS%?(P, Q, K)a a m (P)a b n (QX(K) 

+3igS™f(P, Q, k)a a m (P)a b n (Q)A c (k) 



g 2 



«XJ 2 - («) 2 + 2(A a A a a*a b - A^A'a^) 



+g\A a a a a b a b - A a a b a«<) 
-i 9 -a a m {P) {{K - Q)mUQ)CH(K) + e ahc K m UQ)UK)) 

-i 9 -A a ,{p) ((K - Q)oUQ)Zh(K) + e abc K,UQ)UK)) 

+j («, + 2^S< + A a A a )(e H + C) + j%Ua a m a a m + 2A a a a ) 

+^(& + C? + ^oUe H + &) + igea b cPmci(P)a b m (Q)c c (K) 
+ige abc P cl(P)A b (q)c c (K), (6) 
where the symmetrised trilinear coupling is 

S2 n J(P, Q, K) = l -e abc [(P - Q) t 5 mn + (Q- K) m 5 nl + (K - P)J lm ] . (7) 

The second piece is quadratic and corresponds to the T = and the thermal 
counterterms: 

C CT = \{Z^m 2 + (Z. - l)m 2 )^ + {Z%Z X - 1) 

+ 1 -{Z A - IK (K 2 5 mn - (l - £) K m K^ < 
+ 9 j(Z A Z*Z 2 g - l)$y m + \{Z A - \)AltfA% 




(8) 



The multiplicative and additive renormalisation constants are defined from 
the relations between the renormalised and the bare couplings, as listed next: 



The quadratic approximation to the counterterms is sufficient for our calcu- 
lation, since they contribute only at one-loop. 

The new couplings, emerging after the 4-d (T=0) renormalisation con- 
ditions are imposed, are finite from the point of view of the 4-d ultraviolet 
behavior, but should be considered to be the bare couplings of the effective 
field theory, since they explicitly might depend on the 3-d cut-off. 

3 Non-static fluctuations 

3.1 The potential (local) term of the effective action 

The 1-loop contribution to the potential energy of a homogenous "^-background 
is expressed through the standard sum-integral (with the n = mode ex- 
cluded from the sum) 



The meaning of the primed integration symbol and the characteristic features 
of this integral are detailed in Appendix A. Also the counterterms evaluated 
with the background $ belong to the 1-loop potential. The complete 1-loop 




(9) 




(10) 
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expression of the regularised bare potential is given by 




A A , 2 
In — m 



(for the meaning of the notation J 20 , and others below, see Appendix A). Its 
temperature independent cut-off dependences are cancelled by appropriately 
choosing the coefficients of the renormalisation constants in the "classical" 
counterterm expression: 



V, 



« = \{5m\ + Z^m^l + (2Z„ + Z M ) A$* 



(12) 



The list of two-loop diagrams and the algebraic expressions for each of 
them is the same as in case of the effective potential calculation. They appear 
in several papers [IT, 12" in terms of two standard sum-integrals, in addition 
to the I 4 function defined above. These functions in the present case are 
analytic in the propagator masses mf because the zero-frequency modes are 
excluded from the summation: 



H 4 (m 1 ,m 2 .m 3 ) 



L 4 (mi,m 2 ) 



fpi Ip2 Ip3 S(Pl + P2 + P 3 ) 
1 



f f (PjP 2 ) 2 

jpi jp2 Pt{Pt + m{)Pi{Pl + mi 



(13) 



The latter expression cancels algebraically from the sum of the 2-loop contri- 
butions. The main properties of the functions H4 and L4 appear in Appendix 
A. The sum of the two-loop contributions without the counterterms, dis- 
playing typical three-dimensional linear and logarithmic cut-off dependences 
leads to the following regularised expression (for the values of the constants 
K_, and see Appendix A): 



V* 



(2) 
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~256tt 4 ~~ 512tt 4 ~ 128tt 4 

9/20 \ 



01 



+ 



8tt 2 
% 2 A 



A ; 



2 1 



128tt 4 + 32tt 4 



20 3-K"oq 
4 

27<? 2 A 



" 365g 4 27g 2 A A 2 1 1 
_1024tt 2 + 256tt 2 ~ 32^ J J 



(14) 



From this expression beyond the constants, also terms proportional to m 2 $Q 
are omitted. For the required accuracy m 2 can be related on the tree level 
to the T = Higgs-mass. Since this mass is proportional to A, the terms 
proportional to m 2 are actually 0(g 4 \, g 2 \ 2 , A 3 ). 

One also has to compute the 1-loop non-static counterterm contributions, 
whose sum is simply: 



V, 



9 



CT 



m 2 a (Z m + 2Z gl )h(m a ) 



+ \ (^l{Zxi + Zn) + 5mlj I A (m a ) 
+ ^ Q$o(^ai + Z m ) + 5mj J h{m a ). 



(15) 
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The renormalised potential term of the effective theory can be determined 
once the wave function renormalisation constants are known, by choosing 
expressions for 5m 2 and 5Z X to fulfill some temperature independent renor- 
malisation conditions for the potential energy. So, we need the wave function 
rescalings to 1-loop accuracy first. 



3.2 Kinetic terms of the effective theory 

The effective kinetic terms are extracted from the gradient expansion of ap- 
propriately chosen two-point functions. More closely, they are determined 
by the coefficients of the linear term in the expansion of the 2-point func- 
tions into power series with respect to p 2 . Two kinds of diagrams appear at 
one-loop level. The tadpole-type is momentum independent. The bubble di- 
agrams are listed in Appendix B accompanied by the corresponding analytic 
expressions, expanded to p 2 terms. 

One adds to the corresponding analytic expressions the "classical value" 
of the counterterms. The renormalisation constants Z are fixed by requiring 
unit residue for the propagators: 



25g 2 281g 2 
Za ~ 1 + 48^° " 720^' 

7 - 1 + ^ D ^ 



(16) 



(-D — ln(A/T) — In 2n + ^ E ). The gauge coupling renormalisation is found by 
requiring that the gauge mass term proportional to <3>q, generated radiatively, 
vanish (in other words the coupling in front of (A"$ ) 2 stays at g 2 ): 

Z -1 A3d2 D + 70192 + X (17) 
^- 1 ~96^ + 1440^ + 48^" (17) 

In terms of the renormalised fields the kinetic term of the effective La- 
grangian can be written down by using its invariance under spatial gauge 
transformations in the usual form: 

L km = -F^F« + ^(V^V^, i,j = 1,2,3. (18) 
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3.3 Renormalisation 



Now, we can write after applying ([16]) to the $o - fi e ld that form of the 
regularised potential energy, where the remaining cut-off dependence of 4- 
dimensional character is due exclusively to the mass and scalar self-coupling 
renormalisation: 
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(19) 



The final step is to fix the parameters of the renormalised potential energy 
by enforcing certain renormalisation conditions. We are going to use the sim- 
plest conditions, fixing the second and fourth derivatives of the temperature 
independent part od the potential energy at the origin: 



d?V {T —independent) 



d$i 



—m 



d^V (T —independent) \ 

3*3 



(20) 



One pays for this the price of having more complicated relations to the T = 
physical observables. The connection of the Higgs and of the vector masses 
as well as the vacuum expectation value of the Higgs field to the couplings 
renormalised through the above conditions are given in Appendix C. 

The renormalised potential term of the reduced model is finally given by 
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32tt 4 
(21) 



The last term of ( pT| ) can be split into the sum of a finite and an infinite term 
by introducing a 3d scale ^ into it. It is very remarkable, that the 3d scale 
dependence at this stage has a coeffcient which is of opposite sign relative to 
what one has in the 3-d SU(2) Higgs model. For the O(N) scalar models this 
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has been observed by 0. It cannot be accidental but we did not pursue this 
phenomenon in the present paper. The scale /Z3 should not affect the results 
in the exact solution of the 3d model, but might be tuned at finite order, if 
necessary. 

The finite, /13-independent piece of the two-loop thermal mass can be 
parametrized as m(T) 2 = kig 4 + k 2 g 2 X + k^X 2 and the numerical values of 
the coefficients are 

k! = -0.0390912288 
k 2 = -0.0116685842 

k 3 = 0.0027102886. (22) 
3.4 Non-static nonlocality 

The higher terms of the expansion in the external momentum squared (p 2 ) of 
the bubble diagrams of Appendix B give rise to non-local interactions of the 
static fields. For each field a kernel can be introduced, what we shall denote 
by I H (p 2 ), I G (p 2 ) and I A (p 2 ), respectively. They are defined by subtracting 
from the full expressions of the corresponding bubbles the first two terms of 
their gradient expansions, which were taken into account already by the wave 
function renormalisation and the mass renormalisation locally (see 3.2): 

I z (p 2 ) = I z ( p 2 ) - I z (0) - P 2 I Z '(0) (23) 

(Z = H,G,A). These kernels are used in the calculation of the Higgs- 
potential from the effective 3-d model at the static 1-loop level (see section 
5). Their contribution is of the form of Eq.(p9|) in Appendix E. In place 
of the symbol m in Eq . (|89|) one should use the mass of the corresponding 
static field. When working with 0(g A , Xg, X 2 ) accuracy, we should use the 
approximation (|95D to the p-integral. 

Since I z (0) = 0, the first term of fl9"5| ) does not contribute. This is 
welcome, since this circumstance ensures that no potential linear in $ will 
be produced. Also, we notice that I z (p 2 ) = lf{p 2 ) (for the notations, see 
(pQ"D), therefore the integrands of the last terms are given with help of a single 
function. Nonetheless, as explained in Appendix E, in the second term on 
the right hand side of ( |95D one can use the expansion of I z (p 2 ) with respect 
to the mass-squares of the non-static fields truncated at linear order, while 
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in the third term only the mass-independent piece is to be retained. The 
expression to be used in the first integral we shall denote below by I Z1 , while 
the second by I Z2 . 

One can point out few useful features of the integrals appearing in (|95|), 
allowing the omission of unimportant terms from the full expressions of the 
kernels. We shall discuss these simplifying remarks first, and present only 
the relevant pieces from the expressions of the kernels. 

There are two kinds of diagrams contributing to each I z : 

i) the vertices are independent of the background $o, and are proportional 
to g, 

ii) the vertices are proportional to $o and to g 2 , A. In this case the two- 
point functions actually correspond to certain 4-point vertices, involving two 
constant external $-lines. The corresponding non-locality has been discussed 
for the N- vector model in H. The way we handle them in the present paper 
seems to us more standardizable. 

In case i) the first term of the mass expansion of I ZI is going to contribute 
to the effective Higgs-potential a constant which will be omitted. On the 
other hand kernels of this type when multiplied by m 2 in front of the third 
term of ( j95|) already reach the accuracy limiting our calculation, therefore 
their mass independent part is sufficient for the computation. 

In case ii) the coupling factors in front of each diagram are already 
0(g 4 , etc.), therefore they do not contribute to I Z2 in a calculation with 
the present accuracy, while only their mass independent terms are to be used 
in I Z1 . 

Furthermore, some terms of I z (p 2 ) give rise in the calculation of the ef- 
fective Higgs potential to IR-finite p-integrals fully factorised from the rest 
(this is true in particular for the subtractions from I z (p 2 )). These contribu- 
tions turn out to be proportional to the UV cut-off and their only role is to 
cancel against the 3-dimensional " counterterms" generated in the reduction 
step (see Eq.QZI|)). We need explicitly only the finite contributions from the 
integrals. There is no need to present those parts of the kernels which surely 
do not have any finite contribution. 

With these remarks, we restrict ourselves to presenting only the relevant 
part of the kernels for the three static fields (H,G,A), and also for the static 
ghost. The occurring integrals as a rule were reduced to unit numerators, 
with help of the usual "reduction" procedure [|TT|1 . 
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The Higgs-nonlocality: 
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The magnetic vector nonlocality: 
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The ghost-nonlocality (no ghost contribution arises to the second integral 



of 



2_2„2 



K 4 (K + p) 2 ' 



(27) 
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Two useful remarks can be made concerning the magnetic vector nonlocality: 
In the second equation of ( ffif) sometimes only combinations of terms in 
the integrands prove to be proportional to p 2 . Therefore the corresponding 
weighted integrals should not be performed term by term. 

There are parts of the kernel (|^) which could not be reduced to unit 
numerator, but are proportional to powers of k 2 . We shall see that upon 
calculating their contribution to the Higgs effective potential from the 3-d 
gauge-Higgs effective system, they combine with non-local contributions to 
the potential from the static Ao integration (see the end of subsection 4.1 
below) into formally Lorentz-invariant contributions. 



4 Static electric fluctuations 

4.1 Contribution to the potential term 

The Ao-integration resummed with help of a thermal mass nip yields at 
1-loop: 

V {1 - loop) = \h{m A0 ), (28) 
defined through the integral 

h{m) =T [ ln(k 2 + m 2 ) (29) 

(its value appears in Appendix A). 

At two loops, vacuum diagrams involving one Aq propagator are listed 
with their analytic definitions in Appendix D. In addition also three coun- 
terterm contributions of @ should be taken into account on 1-loop level. 
They can be expressed (up to constant terms) with help of I^uiao) where 
the prime denotes differentiation with respect to m\ : 



V CT ] = ^- J 3( m A0) 



# 2 $ 2 

-m 2 D - Z Al m 2 M H ~r-{.Z A1 + + 2Z gl ) 



(30) 



When only terms up to 0(g A ) are retained and further constants are thrown 
away, one finds the contribution 

(ao) 3Tm 2 n m A o 3o 2 T^A, . , . 

v (ao) = _^_4o + y_j {Zn + 2ZgiY (31) 
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The evaluation of 2-loop diagrams with topology of 'Figure 8' does not 
pose any problem, since their expressions are automatically factorised into 
the product of static and nonstatic (sum-)integrals. The evaluation scheme 
of 'setting sun' diagrams needs, however, some explanation. 

Their general form is as appears in (jB9f). Exploiting the analysis of 
Appendix E we see that the final result of the p-integration depends non- 
analytically on m\ . The non-analytic piece comes from the first term of the 
Taylor expansion of the kernel I with respect to p 2 : 1(0) = I' z (mAo)- in other 
words, this contribution is the product of a static and of a nonstatic integral, 
similarly to the 'Figure 8' topology. In summary, we find that contributions 
linear in m^o come from the thermal counterterm and from the "factorised" 
2-loop expressions. They sum up to 

^(4 AT+ ^f_ m 2). (32 ) 
8tt \ 2?r 2 6 D j K ' 

An optimal choice for the resummation mass (mo) is when it is not re- 
ceiving any finite contribution from higher orders of the perturbation theory. 
This requirement leads to the equality 

m 2 D = h -g 2 T\ (33) 

which means that the two- loop A -contribution at the level of the "local 
approximation" yields exclusively a linearly diverging, non-polynomial con- 
tribution to the potential energy term of the effective 3-d gauge-Higgs model: 

Vt l °° p) = ^4-aoAT. (34) 

This term plays essential role in the consistent 2-loop perturbative solution 
of the effective model (see Section 5). 

Further terms of the expansion of I(p 2 ) are interpreted as higher deriva- 
tive kinetic terms of Aq. Their evaluation is based on (pop. The only dif- 
ference relative to the content of subsection 3.4 is that the A integration is 
performed already at this stage, therefore the contribution from its non-local 
kinetic action modifies now the potential energy of the Higgs field. When 
the terms yielding 0(g 4 , g 2 X, A 2 ) accuracy are retained, it turns out that they 
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actually contribute only to the mass term: 
$t$ J_ 




01 
Q 6 

Q2 + 2Ql + Qiy^_ 



I 



p 2 Q 2 K* 8" J p 2 Q 2 K 2 

f — L_ (V - 2^! + Q 2 ^l _ 
2 7 p 4 Q 2 ^ 2 \ ZP 1 K 2 + Qo K 2 Q 2 Wo J 

99 4 ( Ql{K 2 -Q 2 ) \ 
'~2~J p*Q 2 K* I 

In view of the remark at the end of subsection 3.4 on non-static nonlocali- 
ties, it is not convenient to evaluate explicitly this contribution already at the 
present stage. Its Lorentz non-invariant pieces are going to combine in the 
expression of the final effective Higgs potential with contributions from the 
static magnetic vector fluctuations into simpler Lorentz invariant integrals 
(see section 5). 

4.2 Aq contribution to the gauge-Higgs kinetic terms 

The A -bubble contributes to the self-energy of the magnetic vector and of 
the Higgs (SU(2) singlet) scalar fluctuations. The corresponding analytic 
expressions are given in Appendix B. 

There are two ways to treat these contributions. The first is to apply 
the gradient expansion again, and retain the finite mass correction and the 
wave function rescaling factor from the infinite series. This is the approach 
we followed in case of non-static fluctuations. Then the magnetic vector 
receives only field rescaling correction: 

sz ™ - & < 36 » 

The Higgs particle receives both finite mass correction and field rescaling 
factor from the A -bubble: 

H _ 3g^%T 
~ _ 64^V 
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One should note that the rescaling factors are 0(g) what is just the order 
one has to keep to have the 0(g A ) accurate effective action |TJ|. The different 



behavior of the Higgs and of the Goldstone fields reflects the breakdown of 
the gauge symmetry in presence of $o- 

The Higgs-field mass and field rescaling corrections appear only in a non- 
zero $o background. The $ -dependence of these quantities will be treated 
as a non-fluctuating parametric effect in the course of solving the effective 
gauge-Higgs model. 

In this case the kernel of the one loop integral is not analytic in p 2 , 
therefore no basis can be provided for the gradient expansion. Therefore it is 
important to proceed also the other way, that is to keep the original bubble as 
a whole in form of a nonlocal part of the effective action. Since its coefficient 
is g 4 , in the perturbative solution of the effective model it is sufficient to 
compute its contribution to 1-loop accuracy. (This amounts actually to a 
two-step evaluation of the of the purely static AAH and AAV setting sun 
diagrams.) 

The difference between the two approaches from the point of view of the 
quantitative characterisation of the phase transition will be discussed within 
the perturbative framework in the next Section. 



5 Two-loop effective Higgs potential from the 
3-d effective theory 

In this and the next sections we are going to compare various versions of the 
3-d theory by calculating perturbatively with their help the characterisation 
of the electroweak phase transition. This we achieve by finding the respective 
point of degeneracy of the effective Higgs potential in each approximation. 
The analysis will be performed for g = 2/3 and M w = 80.6GeV, and for 
various choices of the T = Higgs mass. (The scheme of the determination 
of the renormalised parameters is sketched in Appendix C). 
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5.1 Approximations to the effective theory 



We start by summarising the effective Lagrangian obtained from the calcu- 
lations of sections 3 and 4. The parameters were calculated with accuracy 
0{g\\g\X): 



C 



3D 



where 



^^•+^(V l «f) t (V J $)+^m(T) 2 <|.t$+^($t$)2_i- m 



AO~^~£'nonloc~\~£'CTi 

(38) 



m 



AO 



m(Tf 



9i-A. n — 9jA m + g^abcA-i Aj 

d % - ig 3 A^r a , 



9l 



-T + 

6 4/ 



m 2 + T 2 



3 o A a 
IS 9 + 12 +9 



32 



'20 



+</- A ( -ho + 2-^00 

\k 00 



8tt 21 16 



Sg 2 \ 



192tT 



Anonloc 
AO 

3 ' T 



81 157 
16 00 + 7680^ 



(39) 



The dependence of the thermal mass on the 3-d scale is incorrect. The only 
way to ensure the correct scale dependence of the thermal mass is to include 
also into the local approximations the effect of the non-static and and some 
of the static nonlocalities. 

For this we write down their contribution to the effective Higgs potential 
combined from (24), (|25|) (p7[) on one hand, and from (|26"D and (|35|) on 
the other. The contributions from the different nonlocalities are expressed 
through some integrals appearing in Appendix A: 

Higgs+Goldstone: 



1 



if 43 (0,0,0)(-A 2 -^ 4 + 3A^ 



d 



•15 



dm{ 



C(0,0) (™g A + \W 



(40) 
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Magnetic and electric vector: 



2 



^43(0, 0, 0) (±§V + \\g 2 ) + -JL C (0, 0) - f A, 2 ; 



^ 4 rL 4 (0,0) 



dm 2 



(41) 



Ghost: 



2 



^ 43 (o,o,o) + ^ 4 ^ ! c(o,o) 



(42) 



Using the features of the 3 basic integrals (H^,C, L4), listed in Appendix 
A, the final numerical expression for the sum of the local (e.g. ([22])) plus 
non-local 0(g A ) thermal mass terms is the following: 



1$2 

2 



(| , 0.0918423^| + 0.0314655A 3 ^ 



0.0047642A2 



1 In Bl 

16tt 2 111 T 



3 \n 2 81 a 



(43) 



The 2-loop contribution to m 2 (T) could be much influenced by the choice 
of the 3d renormalisation scale. Through the 4-d renormalisation of the mass- 
parameter m? it is sensitive also to the 4-d normalisation scale (see Appendix 
C). 

All fields in ([38]) are renormalised 4-d fields multiplied by y/T. The 3-d 
gauge coupling in the kinetic terms is simply g 3 = gy/T, due to the fact that 
all finite T- independent corrections are included into Z g (cf. eq.flTTD). Also 
the quartic coupling takes the simple form A3 = AT, in view of the particular 
renormalisation condition (f20|). 

The cut-off dependent part of the mass term is hidden in Lqt- It ensures 
the cancellation of the cut-off dependence of the 3-d calculations. Except for 
the non-polynomial term, we are not going to discuss this cancellation. The 
formal correctness of our effective Higgs potential will be verified by checking 
its agreement with results of other 3-d calculations, expressed through the pa- 
rameters m 2 (T),gl, A3, rrijj, tuq, m 2 a . Off course these parameters themselves 



are connected to the 4-d couplings differently. 

The non-polynomial term of the above Lagrangian is simply the result of 
the 1-loop y4 -integration. 
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Terms of the remaining non-local Lagrangian C non i oc are all of the generic 
form 

Cnonioc = -<P(-P)N(P)<P(P) (44) 

(4>(p) refers to general 3-d fields). The kernels M{p) are the result of purely 
static v4 Toops, listed in Appendix B. 



K a 1 (p) = -^g 2 k 



3 ff3rF,2 f t 

k (P + mi )((fc + p)2 + mi ) 



(F + m\ )((k + pf + mi ) 3 (F + m * „)2 



(45) 



In order to see clearly the effect of the non-polynomial and nonlocal terms 
we shall consider the effective potential in three different approximations. 

i) Local, polynomial approximation 

In this approximation non-local effects due to the A -bubbles are ne- 
glected. The parametric "susceptibilities" of the Aj and fields arising 
from the static A -bubble are included into the approximation (the back- 
ground $o appearing in their expressions is a non- fluctuating parameter!). 
The non-polynomial term is expanded to fourth power around m? D . In pres- 
ence of non-zero Higgs background the Lagrangean actually breaks(!) the 
gauge symmetry (some coefficients explicitly depend on the value of the back- 
ground field, cf. ([37])). The approximate form of the effective Lagrangian 
is: 



£lp = ^F^ + \{W^)\W^) + l -Z H {d t i H 

V 3 (&$) + ig 3 e abcPt cl(p)A*(q)c c (k), (46) 



with the local potential 

V 3 (&$) = \(m(T) 2 + 5m 2 )$t$ + X -8m\il + ^^($ f $) 2 , (47) 



where 

" ~~ ' ' 7n : I ■ "A = 

24vrmV ' 64vrm D 



- n 4- \-i x\ _ 9 ^3 

A* - l 1 + 94-Trm -J ' ^ A 
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z " = siS- < 48 > 

(the expressions of m(T) 2 , c/f, A 3 have been given above). We note, that 
also a term Z H (A°:) 2 ^ 2 H /2 is present, but it would contribute to higher orders 
in the couplings, and such terms will not be displayed here. 

In order to apply perturbation theory with canonical kinetic terms, one 
has to rescale both the Ai and fields: 



93 



= g^l\ A i = A i »- 1 ' 2 , l=i H {l + Z H /2). (49) 



Below we continue to use the ■unbarred notation for these quantities! 

After rescalings and the $ -shift in the i? a -gauge the action density is 
written as 



Clp = V 3 ((l-Z H )$l) 

1 A«(k) + clk 2 c a (k) 



+\^{-k) 



(k 2 + ml)Sij - (1 - \)hkj 



+ k H (-k)[k 2 + m 2 h }^ H (k) + l -U-k)[k 2 + ml]Uk) 



+i9sSf c (p, q, k)AUp)A%q)At(k) + ±A*\g*& + 

+f [(A°A) 2 - (A- a-) 2 ] + ^A 2 ^ 

+igze ahc p i c\(p)A a i (q)c c (k) 

-i2 A 1{p)[9h{k - q)i^ a (q)^H(k) + gze abc k^ b {q)^ c {k)\ 

*o^h(?W£ + QhggC) + + Y^^aoA- (50) 

Here the expressions for the different couplings read as 

9h = 9K 1 ~ z h), A gg = A 3 + 5A 

Xhh = (A 3 + SX)(1 - 2Z H ), \ HG = (A3 + 5\){l - Z H ), 
q H HH = (\3 + S\)(l-2Z H ), q HGG = (\ 3 + 5\)(l-Z H ). (51) 



23 



The effective masses of the different fields are 
ml = && (1-Z H ), 

m 2 h = ( m (T) 2 + 5m 2 + 5m 2 H )(l-Z H ) + ^±^<l> 2 (l-2Z H ), 

m 2 g = m(T) 2 + 6m 2 + ^^e (l-Z H ). (52) 

ii) Local, non-polynomial approximation 

The rescaling and the shift of the fields proceeds as before. The main 
difference is that here the order of expanding and shifting the <3>-field in 
the non-polynomial term is changed relative to case i). The form of the 
Lagrangian can be equally written in the form (0). Some expressions and 
constants appearing in it are modified relative to case i): 

V 3 (&$) = ^(m(T) 2 + W)$t$ + ^m^ + ^^($ t $) 2 -^mi , (53) 



5m 2 



ni 



h 



m„ 



m D + iQhK 
3g 2 m A0 



16n 

(m(T) 2 + 5m 2 + 5m 2 H )(l 

m(T) 2 + Sm 2 + ^$ 2 (1 - 
o 



2Z 



H , 



(54) 



Q-HGG — ^3(1 Zh] 



Qhhh = Qhgg ( 1 — Z H 



64:71171 AO ' 

256m A0 ' 



64:7rmAo 



\hg - Xaa(l - Z H) -r 



1 



\hh = A HG (1 - Z H ) + aft » ^ " . (55) 

GV ; 2567rm| 1024vrm^ v ; 
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Hi) Nonlocal, nonpolynomial representation 

Its Lagrangean density coincides with flSSD, and the major issue of our in- 
vestigation is the quantitative comparison of its 0(g A ) solution to the approx- 
imate solutions corresponding to cases i) and ii). The couplings quHH, Ihgg 
and Xhgi ^gg, ^hh agree with their expression derived for case ii) when one 
puts Z H = 0, 5m 2 H = 0. 

5.2 Variations on the effective Higgs potential 

The 1-loop correction to the "classical" value of V3($q) has the unique form 
for all approximations: 



In this expression cut-off dependent terms with polynomial dependence on $ 
are not displayed, their cancellation is taken for granted. Though the form of 
this expression is unique, one should keep in mind the variable meaning of the 
different masses from one approximation to the other. From the point of view 
of consistent cancellation of non-polynomial divergences it is important to 
notice, that the 1-loop linear divergences of the Higgs- and Goldstone-fields 
by the contribution of the shifted non-polynomial term to their respective 
mass produce a non-polynomial cut-off dependent term (— 3g 2 /167r 3 )Am/io 
(in cases ii) and iii) only) to be compared with the induced non-polynomial 
divergence eq. fl34|) . In the first two approximation schemes no further di- 
vergences of this form are produced, therefore their cut-off dependence is 
clearly different from what is generated in the course of the integration. Be- 
low, for cases i) and ii) we assume the presence of the appropriately modified 
counterterms. 

The 2-loop diagrams of a local, 3-d gauge-Higgs system are the same as 
in 4-d. Just the functions appearing in the expressions are defined as three- 
dimensional momentum integrals. In purely algebraic steps one can work out 
the general representation given in terms of the functions I^m 2 ), H 3 (mi, m 2 , m^) 
(see Appendix A), to arrive finally at the following 2-loop contribution: 




l -(Qml + ml + 3ml) + ^m A0 A + V^f n {%). (56) 



v; 



,(2) 
eff 




Q.GGH + 1HHH ^,2 

12 $0 




1287T 2 



(5X G Gm G + 2\ GH m G m H + \ H Hm 2 H ) 
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J^((9-631n3)m 2 - 



128tt 2 V v ' a g 2 

Qg 2 h m a m H 2 t>g\m G m H 3glm 2 H 
- h jm G H g 1 - 2 



3flM . 10 3^m a m G 
- + 12m a m G 



3# 2 m H -m G 2 



G i 



128vr 2 m a 

1 no„2 T 1927T 2 ^ 



1927T 2 

3 



T 



64vr 2 



2 2 l ''''a 

yh H 2m a 



■a + ra g ^I^h ln m H (2m a + m H ) 
4m 2 



2m 2 



+ m H 
, + m 
TfiQ + m H 



(m a + m H ) 



m 2 G ) 2 m a + m G + m H 



-g 2 h (m 2 H + m 2 G -^)\n 



m a + m G + m H 



T 



- 2m G ) In 



m a + 2m G 



T 



gjmj , m a + m H 



T 



+2g 2 h m 2 a ln 



2m a + m H 
T 



— 11m 2 ln — ^ 



T 



(57) 



The final step of composing the full effective potential corresponds to picking 
up the contributions of the purely static VAA, and HAA diagrams: 



V. 



(2nonloc) 



eff 



6g 2 L (m 2 A0 - -m 2 ) + —^{m 2 M + 2m A om a 



M 

327T 2 



T7rl Am 



,40 



3<7 3 2 



— m 2 ln 

2 a 



3 [(4mi - m 2 a ) ln 
2m A0 + ^i/ n 



2m A o + m a 



(58) 



The cancellation of the linear divergence, nonpolynomial in $q can be seen 
explicitly. 

The 2-loop corrections to the effective potential, given by Eqs. ([57D and 
(|58|) can be compared to the results of the direct 4-d calculation of [|l^, 
Those calculations were done with different regularisation and renormalisa- 
tion scheme, therefore in the expression of the potential terms reflecting this 
difference are seen. Still, using the leading expressions for our couplings 
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9hi Igghi Qhhh, Aggj ^gh and Xhh, all logarithmic terms expressed through 
the propagator masses agree. Also our polynomial terms, except those pro- 
portional to the regularisation dependent coefficient L Q , have their exact 
counterpart in the 4-d expression. The T-independent contribution appear- 
ing in the 4d result has no corresponding term in our case, an obvious effect of 
the difference of the renormalisation schemes. The polynomial terms explic- 
itly depending on jx and proportional to the constants ci, c 2 characteristic for 
the dimensional regularisation, can be compared to our terms proportional 
to Lq. Choosing in both potentials n = T, one finds in the 4-d expression 
for the coefficients of the terms g 2 m 2 v , g 2 (m 2 H + 3m G ) and A 2 $o, the values 
2.441 x 1(T 2 , 4.012 x 10~ 3 , -2.239 x 10~ 4 , respectively. The corresponding co- 
efficients from our expression are: 2.765 x 1CT 2 , 5.027 x 10~ 3 , —5.374 x 10" 4 . 
Therefore the only origin of considerable deviations in physical quantities 
could be the effect of the reduction on various couplings. Our numerical ex- 
perience was that the 0(g A )T 2 correction in m 2 (T) accounts for essentially 
all differences. 



6 Phase transition characteristics 

In this section we describe and discuss the phase transition in the three per- 
turbative approximations introduced in the previous section. Our analysis 
will cover the Higgs mass range 30 to 120 GeV. Our main interest lies in find- 
ing the percentual variation in the following physical quantities: the critical 
temperature (T c ), the Higgs discontinuity (<& c ), the surface tension (a) and 
the latent heat (L c ). The amount of variation from one approximation to 
the other will give an idea of the theoretical error coming from the reduction 
step. 

The first step of the calculation consists to choose values for g and A. For 
example one can fix g = 2/3 and tune A guided by the tree-level Higgs-to- 
gauge boson mass ratio. All dimensional quantities are scaled by appropriate 
powers of \m\ (what practically amounts to set \m\ = 1 in the expressions of 
the Higgs effective potential). 

Next, one finds from the degeneracy condition of the effective potential 
the ratio T c /\m\. Here we have to discuss a phenomenon already noticed by 
[ |l!| . It has been observed that for Higgs mass values ra# > lOOGeV when the 



temperature is lowered one reaches the barrier temperature before the dege- 
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naracy of the minima would occur. The phenomenon was traced back to the 
action of the term ~ g 2 mAo{ , mH ^^mc) in the effective potential. In our non- 
polynomial approximation the y4 -integration contributes a negative term to 
m 2 (T) (i.e. 5m 2 ), which acts even stronger, and the same phenomenon is gen- 
erated, when using // 3 = T, already for m H < 70GeV. However, by choosing 
somewhat more exotic value for the normalisation scale (fi^ = T/17), we 
could follow the transition up to m H = 120GeV. We have used this value in 
the whole temperature range. It has been checked for m# = 35GeV that the 
variation of /i3 leads in the physical quantities to negligible changes. 

Finally, the relations of Appendix C allow to express with 0(g 4 ) accuracy 
the ratios Mw/\m\ and Mn/\m\ with help of g, A and T c /\m\ (the dependence 
on the latter appears through our choice of the 4-d normalisation point). 
The pole mass M# resulting from this calculation appears in column 2 of 
Table 1, where M w = 80.6GeV sets the scale in physical units. Our 4- 
d renormalisation scheme leads to somewhat smaller mass shifts, than the 
scheme used in fll~5|. 



In order to present physically meaningful plots one has to eliminate from 
these relations \m\ and scale everything by an appropriate physical quantity. 
In Fig. 1 (coulumn 3 of Table 1) we present T c /Mh as a function of Mh/Mw 
in the non-polynomial non-local (NNA) approximation (upper case notation 
for the masses always refer to T=0 pole mass). Internally, our different 
approximations do not affect the critical temperature, they all agree within 
better than 1%, as can be seen in Table 2. The curve in Fig.l is, however, 



systematically below the data appearing in Table 2 of [TE\. 15% deviation 



is observed for m# =35GeV, which gradually decreases to 8% for m H 



90GeV. If one wishes to compare to the 2-loop 4d calculations of |L2], |14 
one has to use their coupling point: g = 0.63722, mw = 80.22. Again we 
find that our T c /rriH values are 10% lower than those appearing in Fig. 4 of 
[0. Our 0(g 4 )T 2 correction to m 2 {T) is about 10% of the 1-loop value 
(and is about 9 times larger than what was found in |TJ|), therefore one 
qualitatively understands that the barrier temperature is brought down in 
about the same proportion. At least in the region Mw — Mh the transition 
is already so weak that the transition temperature should agree very well 
with the temperatures limiting the metastability range. 

In Fig. 2 (column 4 of Table 1) the order parameter discontinuity is 
displayed in proportion to T c . Here the agreement is extremely good in the 



whole 35-90 GeV range both in comparison to [T5|, [J and |TJ]]. Also the 
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variance of our different approximations (see Table 2) is minimal. 

The most interesting is the case of the surface tension, which is shown for 
all three approximations in Fig. 4 in (GeV) A units (the dimensionless ratio 
appears in column 5 of Table 1). We did not observe any strengthening ten- 
dency for larger Higgs masses, in contrast to |12[]. This systematic difference 
seems to be correlated with the extended range where we find phase transi- 
tion. It leads to o c values in the range Mh = 50 — 80GeV which perfectly 
agree with the perturbative values quoted by [f|. The dispersion between 
the values of our different approximations is much larger in the high mass 
range (Table 2) than for other quantities, reflecting the increased sensitivity 
to non-local and non-polynomial contributions. 

The situation is just a little less satisfactory in case of the latent heat 
(Fig. 3 and column 6 of Table 1). Our approximations are 10-15% above 
L c /T^ curve of Fig. 6 of and also of the perturbative values in Table 6 of 



7 Conclusions 

In this paper we have attempted to discuss in great detail the reduction strat- 
egy, allowing non-renormalisable and non-local effective interactions. We 
have made explicit the effect of various approximations on the phase tran- 
sition characteristics of the SU(2) Higgs model. Our investigation remained 
within the framework of the perturbation theory. By comparing with other 



perturbative studies |L2j, [14], [L5| we have pointed out that the T C /M H ratio 
is quite sensitive to the choice of the 3-d renormalisation scale. The order 
parameter discontinuity and the the surface tension were found only moder- 
ately sensitive, while the latent heat in the MS scheme seems to drop faster 
with increasing Higgs mass. The minimum of the surface tension when Mh 
is varied has disappeared. One might wonder to what extent is the strength- 
ening effect observed in the 4-d perturbative treatment with increasing Higgs 
mass physical at all. 

The local polynomial and the local non-polynomial approximations start 
to show important (larger than 5%) deviations from the nonlocal, nonpolyno- 
mial version of the effective model only for Mh > 80GeV (also below 30GeV). 
Since in these regions the perturbation theory becomes anyhow unreliable, 
we can say that the application of the reduction to the description of the 
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electroweak phase transition in the relevant Higgs mass range (Mh ~ A%) 
could be as accurate as 5%. This can be true for dimensionless ratios, but 
not for T c /Mh- 

The present assessment of the accuracy is not universal, though the 
structure of the analysis is quite well fitting other field theoretic models, 
too. For instance, our methods could be applied to extended versions of 



the electroweak theory [17], [18|, [19| which are accessible to quantitative non- 



perturbative investigations only in 3-d reduced form. 

One knows that non-perturbative studies of the SU(2) Higgs transition 
led to lower T c /Mh ratio, than the perturbative prediction [p0[] , and at least 
for Mh = 80GeV the surface tension is much lower than it was thought on 
the basis of the strenghtening effect JH]. The results from the PPI-reduction 
seems to push the perturbative phase transition characteristics towards this 
direction. The expressions of the mass parameter of the effective theory (PT|) 
and the T=0 pole masses (Appendix C) provide the necessary background 
for the analysis of 3-d non-perturbative investigations following the strategy 
advocated by [16|, [J. This will be the subject of a forthcoming publication 



A Important integrals 

All integrals have been performed with 3-d momentum cut-off A. The capital 
subscript denotes finite temperature sum-integral, the bold lower case index 
refers to 3-d momentum space integration. The prime on the sum-integral 
means summation over the non-static modes only. 
1. 

h{m) = I + hm 2 + I 2 m 4 + ... 

A 2 AT T 2 

1 ~ 8^ 2 ~27 2 + 12 

= ii 2 A 2 + J n AT + I W T 2 

1 A 1 , , , . 

ln ^ + T7rr( 1 + ln ( 27r )-7s) 



167T 2 T 167T 2 

A 
T 



I 2ln \n^ + I 20 . (59) 
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2. 



h(m) = 7^ Am 2 _ 7~ m3 



2J 1 Am 2 + 2J m 3 . (60) 



3. 



H 4 ( mi ,m 2 ,m 3 ) = K + K 1 m ' +m 3 ' + m ' +... 

5 A 

K = 0.0001041333A 2 - 0.0029850437AT + rT 2 ln- 

32tt 2 T 

-0.0152887686T 2 
= K 02 A 2 + K 01 AT + K 0ln \n^ + K 00 T 2 

K x = -^^(ln^) 2 + 0.001087871 In ^ + const. 

= K 10 (\n ^) 2 + K llog In ^ + const.. (61) 



1 A 2 

# 3 (mi,m2,m 3 ) = -— In ■ ■ + L , 

62ir z {mi + m 2 + m 3 ) z 

-3 



6.70322 x (62) 



H. 



l:! ~ 7pi 7p 2 7 P (P 2 + m 2 ) (P 2 + m 2 2 ) (p 2 + m 2 ) 

a.AT - ^T 2 In A + fl2 T 2 + |g(ln £ - 16^) 

+asm 3 _ 768^2 ( m ? + m 2) + ---, ( 63 ) 



ai = 0.000995014555, a 2 = 0.0074744795, a 3 = 0.0000659642. (64) 
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6. 



L 4 (m 1 ,m 2 ) 



7. 



C(m 1 ,m 2 ) 



(P1P2 



pi Jp2 P?(P? + ml) Pi {Pi + ml) 
T 2 



Lao + 
AT 



768tt 2 
A 



1287T 4 



ln--(l + ln(27r)- 7B )) 



(m 2 + m 2 ). (65) 



1 



1 



T ' 5 A 7 q ( W 2 + k 2 + m? ) (u; 2 + q 2 + ml) 
= C + {m\ + m\)C x + 



(66) 



^ = -l^AT(log^-l-M2.) +7a )-^. (67) 

B Momentum dependent (bubble) contribu- 
tions to 2-point functions 

In this Appendix we give a list of the diagrams generating higher momentum- 
dependent vertices of the reduced theory. Each bubble contains two internal 
lines, which we use for "naming" the diagrams. (The abbreviations are: V 
- for non-static magnetic vector potential, H - for non-static Higgs, G - for 
non-static Goldstone, C - for non-static ghost and A - for static electric 
potential). Some diagrams are proportional also to $q, so they actually cor- 
respond to 4-point functions. The 2-point approach is just a convenient way 
to study their momentum dependence. All 2-point functions depend on the 
spatial momentum K = (0, k). The first few terms of the gradient expansion 
of each diagram will be given up to that term which should be considered in 
the evaluation of local mass and wave function renormalisations. These terms 
are subtracted from the full expression, when one turns to the investigation 
of non-local effects (see Appendix E). 
Higgs 2-point function 
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^VG = 

8 



-9 2 k 2 



(KPf 



1 





o P 2 ' 




p~ 2 . 



P 2 ' + ml(p + k) 2 + m 2 G 
1 1 



P 2 + m 2 a P 2 + m 2 G 



(68) 



^-$ 2 VV = 

32 



1 



8 Jp P 2 + m 2 a (P + k) 2 + w 2 
9g*$l f 1 



(P(P + fc)) ; 
P 2 (P + k) 2 



8 Jp (P 2 + ml) 



2\2 ' 



(69) 



-$ 2 HH = 

72 



1 



1 



A 2 $ 2 

2 Jp P 2 + m 2 H (P + k) 2 + m 2 H 
A 2 $ 2 , i 



2 Jp(P 2 + m 2 H ) 2 ' 
A 2 $ 2 



6 

A 2 $ 2 



JpP 2 + 

2 . 



(P + k) 2 + 



^$ 2 AA = 

32 



6 Jp (P 2 + m 2 G ) 2 ' 
3<? 4 $ 2 r 1 1 



(70) 



(71) 
(72) 



8 J P p 2 + m 2 M (p + k) 2 + m 2 A0 ' 

This last integral is not analytic in k 2 , therefore its gradient expansion is a 
dangerous step. 



Goldstone 2-point function 



-VG 



= - 2g 2 J p (k 2 - 
~ -2g 2 k 2 j 



(KP) 2 1 1 

P 2 P 2 + m 2 (P + A;) 2 + m 2 G 



3P 2 



1 



1 



P 2 + m 2 P 2 + 



m 



2 ' 
G 



(73) 
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-9 2 k 2 



{KPf 1 1 

P 2 ' P 2 + ml(P + k) 2 + m 2 H 



p" 
3P 2 



P 2 + m 2 P 2 + m# ' 



(74) 



A 2 $ 



A 2 $ 2 



L 



9 Jp P 2 + (P + k) 2 + m 2 , 
1 1 



9 Jp P 2 + m 2 H P 2 + m 2 G 



(75) 



Magnetic vector 2-point function 

Here we give the contributions to the polarisation tensor Mij. When the 
static Ai legs are contracted, its contribution to the effective Higgs potential 
is of the form (see Appendix E): 



J p 2 a J 



(M,(p)-^(o))n 



(76) 



where ITjj = 5ij —piPj/p 2 - The following bubbles contribute: 



^VV = 

2 



2% p 2 



(K 2 + m 2 )(Q 2 + m 2 

+TO . ( 3 +2 £ + «) 

A; 2 



K 2 



- % 2 P 2 n 



A; 2 



7 



+ - 



(if 2 + m 2 ) 2 

2 r k 2 r 

3 7 iT 2 (ir 2 + m 2 ) 2 ~ 7 (K 2 + m 2 ) 3 

4 /• ^ ^ fc4 

157 K 2 (K 2 + m 2 ) 3 



(K 2 + m 2 ) 5 



5 7 (fT 2 + m 2 ) 4 
1 r A; 4 



La* 



15 J K A (K 2 + m 2 a ) 2 



(K 2 + m 2 ,)(g 2 + m 2 ; ) 
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9' 



GG 



-^-$ 2 VH 
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^CC 



^AA 



(Q 2 + m 2 H )(Q 2 + m 2 G ) 



i 2tt 



/ 



(g 2 + m^)(g 2 + m^) 2 
inLinL 



' 4 



5 Q 2 + m 2 G 



g2 I (K 2 + m 2 G )(Q 2 + m 2 G ) 



(K 2 + 



m 
A' 2 



h) 2 



[K 2 + m 2 G f 
1 



5K 2 + m 2 



4 ^O 11 *? 



(K 2 + m 2 a )(Q 2 + m 2 H ) 
1 



KILKIL, 



K 2 

{K 2 + m 2 H ){K 2 + m 2 ) ^ 3 K 2 J 



if 2 Q s 



A: 2 

3^ J (k 2 + m 2 M ) 2 



1 - 



4 

+ -■ 



.2 7,2 



k 2 + 5 (A; 2 + m 2 A0 ) 2 



(77) 



(iHIi = kj(5ij -PiPj/p 2 )). 
Ghost 2-point function 



2 



■CY = 2g 2 J K [p 2 - 



{PK) 2 \ 1 1 
K 2 ) K 2 + m\ (P + Kf 



(78) 



C Couplings in terms of T = physical quan- 
tities 

At tree level the renormalised couplings and the vacuum expectation value 
of the Higgs-field (v) determine the masses of the field quanta: 



m 2 H = m 2 + f -w 2 , 
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m 2 w = 9 iv\ (79) 

Physical masses are approximated by finding the roots at p 2 = — m 2 of the 
corresponding 2-point functions. For achieving the accuracy 0(g 4 , ..) we need 
also the vacuum expectation value v corrected by 1-loop fluctuations. For 
this we calculate the minimum of the T = renormalised effective Higgs 
potential to 1-loop from the condition: 



dV f 



dv 2 




4AJ- 



20 ) 



2 




(80) 



For the accuracy requested it is sufficient to use in the 1-loop part the tree 
level expressions (|79|). Then the iterative solution is easily found: 

2 6m 2 . 27g\ T 

1 r 27g\ ,3g 2 , m 2 411 , 2m 2 N 

ln (l7r(-^))+4Aln(---)]). (81) 



64vr 2L 2A v 2A v T 2 " v T 2 

For the Higgs and W-masses one makes use of their 2-point functions. 
Here again we use \x = T and neglect the running of the couplings to the 
scale where they are usually measured (say at m w ). (This should be used 
only in the tree level mass expressions, if necessary.) 

The 4-d integrals were evaluated with cut-off regularisation. We have 
reduced all the 1-loop integrals to 1-variable integrals with help of Feynman 
parametrisation. For a set of numerical values of g, A one finds Mjj and Myy, 
the 1-loop accurate masses in proportion to the renormalised mass parameter 
— m 2 to be still functions of the logarithms of m 2 H /T 2 and of m^/T 2 . In the 
arguments of the logarithms the tree level mass expressions can be used. 
This leads to the unique dependence on —m 2 /T 2 . The numerical value of 
the latter for fixed values of g, A is extracted from the degeneracy condition 
imposed on the minima of the effective potential (see section 6). Having the 
value of ln(— m 2 /T 2 ), the integrations over the Feynman parameter can be 
done. 

Below we give the expressions of the 1-loop accurate Higgs and W-masses 
before the 1-dimensional integrals are performed. 
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The 1-loop accurate expression of the Higgs-mass is 

M 2 H = m 2 H + Y,{-m 2 H ). (82) 
The self-energy of the Higgs field is given by 



2 3g 2 m 2 W /n ,, 2 l,^ T 2 *m 2 T 2 
T,(—m H ) = „ A - (9 — 3h + -h ) In — 5 — r In ■ 

v HJ 64tt 2 v 2 ; m 2 w 8tt 2 tt% 



i 2 ^ 87r 2 m 2 

,2 on Q„2_2 / n U 2 l4 



| | 21 

1 + /i 2 (- - 12x + 12a; 2 ) + h A x{ — + — x - 20x 2 + 10x 3 ) 
2 2 2 



+Z CT (-m 2 H ), (83) 
where /i = m H /m w and the one- variable integral Qh is defined as 

Q H (/> 2 ) = C dx 
Jo 

x ln(l - h 2 x{\ - x)). (84) 
The contribution coming from the 1-loop counterterms is 

_ mi, , . w 9o 2 81o 4 . g 2 m? f 

S CT = — ^-(1 + In 2tt - 7b) (— - 4A 2_) - - f . 85 

G 16tt 2V v ; ' A 2 8A ; 32tt 2 v ; 

We give the W-boson polarisation function at p 2 = —rriyy in an even 
simpler form: 

„ , o , g 2 m 2 /59 + 3/i 2 T 2 3h 2 \ 



' 32tt 2 
Am 2 ^ 
24tt 2 



tte[-l + (/i 2 - 2)x + x 2 ] ln[(l - a;) 2 + /i 2 a;] 

ao 



+ ^ (Li^y ( 86) 

°"- 2 { 2 m 2 H ) K ' 



D Mixed 2-loop contribution to the potential 
energy density of the effective theory 



'Figure 8' diagrams and thermal counterterm: 

^AV = 6g 2 [ - 1 2 / ^ 1 - 

2 J k 2 + m 2 A0 J Q 2 + m 2 a 
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/ k 2 + m 2 / 



V Q 2 (Q 2 + ro*)' 



(AG + AH) 



■\m 2 D A 



3_l 

8 
3 



I k 2 + m 2 A0 { 3 J Q 2 + m 2 G + I Q 2 + m 2 H ) ' 



= —-m 



D 



k 2 + m 2 A0 



(87) 



Here the integration signs refer either to 3-d momentum space integration 
over k, or to finite T sum- integrals over Q. 

Mixed setting-sun diagrams: 



f AW 



(p 2 + m 2 A0 )(Q 2 + m 2 a )(K 2 + m 2 a ) 



2-AHG 



4rAGG 



-^$ 2 AVH 



2p 2 



= ~2 9 



~*9 



K 2 



K 2 Q 2 

Ql 



ij> 2 + m 2 A0 )(Q 2 + m 2 H )(K 2 + m 2 G ) ' 
Ql 



(p 2 + m 2 A0 )(Q 2 + m 2 G )(K 2 + m 2 G y 



4^2 



^ACC =3g 2 J 



( P 2 + m 2 A0 )(Qi + mZ)(K* + m 2 H ) 

Ql 



( P 2 + m 2 A0 )Q 2 K 2 



(88) 



In these integrals the static momentum p and the non-static momenta Q and 
K are related via p + K + Q = 0. 



E Dependence of Feynman-integrals involv- 
ing a static propagator on its mass 

We analyze in Sections 4 and 5 integrals which are of the general form 



I(p 2 



P p z + m l 
38 



(89) 



where I(p 2 ) is the result of non-static sum-integrals. This function has a 
well-defined Taylor-series in p 2 , therefore one can form the following set of 
functions: 

W) = /(P 2 )-E^/ W (0), h = h (90) 

with /w denoting the i-th derivative with respect to p 2 . The dependence of 
/(p 2 ) on the non-static mass-squares is also analytic, the coefficients of the 
corresponding series expansions are IR-finite. 

Below we show that the part of the above integral depending non-analytically 
on m 2 can be constructed explicitly. Namely, we prove, that 

I{p 2 ) m 2 2 ~ 2 f I s (p 2 ) 

77^ = "S /(p = ~ m ) + S ( " m ] i^-y (91) 

The first term on the right hand side is the explicit non-analytic piece, which 
can be interpreted as a factorised product of a static and of a non-static 
integral. The second term is IR-convergent. 
The proof is based on the following identity: 

W) h + i{p 2 ) 1 2 I s+ i(p 2 ) ( . 

^ „l „2/„2 i ™2\ " L „2fs4-2U„2 i ™2\ 



p2(s+V) ^p2 _|_ npp2\ p2(s+2) p2^p2 _|_ yy)2^ p2(s+2)/p2 _|_ ^2 

This identity is applied repeatedly starting with the second term on the right 
hand side of the simple identity: 

Hp-) _'(f)_ ml Jtf) 



p 2 + m 2 p 2 p 2 (p 2 + m 2 ) 

Eventually, one arrives at the following representation of the integrand of 

Hp 2 ) 2v W) _ avv_ 2 , s /^(o)^^_ /q^\ 

p 2 + m 2 - 2LA 771 ) p 2(s+l) m ^A m ) s , p 2 (p 2 + m 2)- ^ 

The second term on the right hand side can be evaluated explicitly, when 
this equality is integrated with respect to p. Also the s-sum in it can be 
done, leading to (pi]). 

In our applications the integral we are interested in is a "mixed" static- 
non-static "setting sun". From its two vertices one has at least a factor g 2 , 
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therefore in an 0(g A ) accurate calculation one is satisfied with keeping the 
first few terms of the expansion ( [94]) (m 2 is of the order of g 2 ): 

On the right hand side in the second term it is sufficient for the present 
accuracy to use the expansion of I(p 2 ) to linear order with respect to the 
non-static mass squares. In the third integral for the same purpose one 
simply sets these masses equal to zero. 
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Table l:Critical quantities from the non-polynomial, nonlocal approxi- 
mation (NNA) 

Table 2:Percentual variation of physical quantities relative to the NNA 
approximation. The symbols 5p and 5jy refer to the (LP A — NNA) /LP A 
and (LNA - NNA) jLN A differences, respectively. 
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Figure Captions 

Figure 1:T C /M H ratio as a function of M H /M W in NNA 

Figure 2:$ C /T C ratio as a function of M H /M W in NNA 

Figure 3:L C /T C 4 ratio as a function of M H /M W in NNA 

Figure 4:The surface tension (cr c ) in (GeV) 3 units as a function of 
M H /M W in various approximations, a: NNA, b: LPA, c: LNA 
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Table 1: Critical quantities from the non-polynomial, nonlocal approximation 
(NNA) 
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Table 2: Percentual variation of physical quantities relative to the NNA 
approximation. The symbols 5 P and 5 N refer to the (LP A — NNA) /LP A 
and (LNA - NNA) jLNA differences, respectively. 
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Figure 1: T c /M# ratio as a function of Mn/M\y in NNA 
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Figure 2: $ C /T C ratio as a function of M H /M W in NNA 
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Figure 3: L c /T c 4 ratio as a function of M H /M W in NNA 
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Figure 4: The surface tension (<r c ) in (GeV) 3 units as a function of M H /M W 
in various approximations, a: NNA, b: LPA, c: LNA 
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